CRM-2931 (2003) 



Second and Third Order Observables of the Two— Matrix Model 



M. Bertola tt2 

t Department of Mathematics and Statistics, Concordia University 
7141 Sherbrooke W., Montreal, Quebec, Canada H^B 1R6 

* Centre de recherches mathematiques, Universite de Montreal 
C. P. 6128, succ. centre ville, Montreal, Quebec, Canada H3C 3J7 



Abstract 

In this paper we complement our recent result on the explicit formula for the planar limit of 
the free energy of the two-matrix model by computing the second and third order observables 
of the model in terms of canonical structures of the underlying genus g spectral curve. In 
particular we provide explicit formulas for any three-loop correlator of the model. Some 
explicit examples are worked out. 

1 Introduction 

The present paper is the natural complement of our recent publication [3] in which an explicit residue formula for 
the free-energy in the planar limit (and for arbitrary genus spectral curves) of the Hcrmitian two-matrix model 
was derived. In point of fact the present results, while logically consequent that paper, can be (and in fact are) 
derived without using the explicit formula in 5 . Given the well-known relation of the free-energy of this model 
with the tau-function of the dispersionless Toda hierarchy, we could as well rephrase the entire contents of this 
(and the previous) paper in the language of integrable hierarchies. 

Let us briefly recall the historical and logical setting of the problem. The object under scrutiny is the random 
2-matrix model [SUE]! which has many applications to solid state physics (23] (e.g., conduction in mesoscopic 
devices, quantum chaos and, lately, crystal growthp2]), particle physics [SJ, 2cZ-quantum gravity and string theory 
|13l 1141 13). The model consists of two Hermitian matrices M\,M 2 of size N x N with a probability distribution 
given by the formula 



d/i(Mi,M 2 ) = ^-dAfidAf 2 exp 



(Fi(Mi) + V 2 {M 2 ) -M 1 M 2 ) 



K=\ J=l 



where the functions V% (called potentials) are polynomials of degree di + 1 for simplicity (and definiteness), but could 
be formally extended to power series. The partition function Zn is known to be a r-function for the KP hierarchy 
in each set of deformation parameters (coefficients of V\ or V 2 ) and to provide solutions of the two-Toda hierarchy 
031 IU H] • This model with polynomial potentials has been investigated in the series of paper |SJ |U] where a 
duality of spectral curves and differential systems for the relevant biorthogonal polynomials has been unveiled 
and analyzed in the case of polynomial potentials. In |10| the mixed correlation functions of the model (traces of 
powers of the two non-commuting matrices) have been reduced to a formal Fredholm-like determinant without any 
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assumption on the nature of the potentials and using the recursion coefficients for the biorthogonal polynomials. 
We recall that the biorthogonal polynomials are two sequences of monic polynomials ( 7 and references therein) 



7r n {x) = x n + ■ ■ ■ , a n (y) = y n + ■ ■ • , n = 0,l,... (1-2) 
that are "orthogonal" (better say "dual") w.r.t. to the coupled measure on the product space 

[dxdy ir n (x)a m {y)e-^ v ^+ v ^-^ = h n 6 mn , ft^OVneN (1-3) 



where V\(x) and V-i^y) are the functions in the measure (jl-ip . It is convenient to introduce the associated quasipoly- 
nomials defined by the formulas 

1> n {x) := -^=-K n - 1 {x)e-^ v ^ (1-4) 

<f> n (y) := -^=a n ^(y)e'r^(y) . (1-5) 
V n-n-i 

In terms of these two sequences of quasipolynomials the multiplications by x and y respectively are represented by 
semiinfinite square matrices Q — [QijjijeN* and P = [Pij]i,jeN* according to the formulae 

Xlp n (x) = ^ Qn,mlpm(x) J y(j>n{v) = ^ Pn.n(f>m(y) 

m m 
Qn,m = = P n ,n, if » > TJ1 + 1. (1-6) 

The matrices P and Q satisfy the "string equation" 

[P,Q] = fti • (1-7) 

We refer for further details to [7| |S1 E| where these models are studied especially in the case of polynomial 
potentials. We also point out that the model can easily be generalized to accommodate contours of integration 
other than the real axes and arbitrary (possibly complex) potentials Cj ■ 

The partition function is believed to have a large N expansion (with H = jj and t fixed) according to the 
formula 

-^lnZ JV =^ = ^) + i_^ 1 > + ... . (1-8) 

This expansion in powers of iV -2 has been repeatedly advocated for the 2-matrix model on the basis of physical 
arguments |141 1171 118j and has been rigorously proven in the one- matrix model |21j . In the two-matrix model 
this expansion is believed to generate 2-dimensional statistical models of surfaces triangulated with ribbon-graphs 
|14U13ll2"5] . where the powers of TV -1 are the Euler characteristics of the surfaces being tessellated. From this point 
of view the term J 7 ^ corresponds to a genus tessellation and the next to a genus one tessellation. Remarkably, 
an algorithm for the computation of the subleading terms in the h = expansion is known and also a closed 
expression of the genus 1 correction J 7 ^ |17l I18j . 

The object of [5] was the leading term of the free energy, J 7 ' -*. It is the generating function (to leading order) 
of the expectations of the powers of the two matrices in the model 

(TtM^) = KduxjM + 0(N- 2 ) , (TrM 2 J ) = Jd VJ F^ + 0(N- 2 ) . (1-9) 

It is at the same time the logarithm of the r-function of the dispersionless Toda hierarchy (in the one-cut case) 
|37II38| : this hierarchy falls in a broad sense in the same sort of integrable dispersionless hierarchies called "universal 
Whitham hierarchy" (i.e. dispersionless KdV and KP) and studied in [3T], where formulas for the corresponding 
t- function are provided which are quite similar to that in pj]. However the situation under inspection in this paper 
(and in pj) involves a hierarchy which does not quite fall in the very broad class of ^Tj; there the moduli space 
under consideration was the set of curves of genus g with N marked points oofe, k = X, . . . , N and local parameters 



2 



around the punctures, together with a fixed meromorphic differential of the second kind d£ with poles of degree 
nu + 1 at oofc. The (local) dimension of the space of these data is 3g — 2 + J2k=i( Uk + ■*-)■ 

The object of interest of the present paper is -among others- the explicit computation of the double and triple 
correlators 

TiM l a TiM j b ) , (TrM i a TrM i b TrM fe c \ , k = 1, 2, a, b, c € N. (1-10) 

/ conn \ / conn 

or -which is the same their generating functions 

1 1 \ S 5 

Tr — Tr > =- — — - — — T 



xi - Mi x 2 ~ Ah I conn 5V 1 ( Xl ) SV t (x 2 ) ' 

1 1 \ s s 

Tr — Tr — > = — , N „„ , 



x-M x y-M 2 / conn SV 1 (x)8V 2 ( y y 

Tr 1 Tr 1 \ = -J—-J—T , (1-11) 

Vi - M 2 2/2 - M 2 / con „ dTMyi) oV2(y 2 ) 



1 1 1 \ <S <5 5 

Tr Tr Tr ) = - — — - — — - — — T 



xi- Mi x 2 -Mi x 3 -Mi/ conn SVi(xi) 8Vi(x 2 ) 5Vi(x 3 ) 

1 1 1 \ 6 S S 

Tr — Tr Tr > = - — — - — — , S F 



xi -Mi sa -Mi y-M 2 / conn SVi{xi) SVi(x 2 ) SV 2 (y) 

1 1 1 \ 8 8 6 

Tr Tr Tr > = - — — - — — - — — T 



yi-M 2 y 2 -M 2 x - M x I conn 8V 2 (yi) SV 2 (y 2 ) 8V x (x) 

Tr - Tr - Tr - \ =— - — T , (1-12) 

2/i - M 2 2/2 - M 2 2/3 " M 2 1 cann 8V 2 (yi) 8V 2 (y 2 ) 8V 2 (y 3 ) 

where the "puncture" operators are defined (a bit formally) as 

8Vi{x) ^ 8V 2 {y) 

The meaning of the puncture operators when the potentials are polynomials is given by the understanding that 
each derivative -say- w.r.t. uk gives a TrMi^ contribution inside the expectation value. The two-loops correlators 
(|1-11|) are well studied and our results here are not new at least in the one-cut case ^2], but are needed for 
computing the three-loop correlator (|1-12H . The latter have not been computed for this model; the similar problem 
for the normal matrix model has been addressed in jl^l in a setting that is a subsetting of our, primarily due to the 
fact their case corresponds to a spectral curve of genus plus certain reality conditions. Indeed the planar limit 
of the free energy !F^> (we will drop the superscripts from now on) can be defined for any spectral curve E s of 
genus g as a closed differential on a certain moduli space and this differential can be integrated explicitly ([5] and 
references therein). In the literature on matrix models the genus of this spectral curve is related to the number 
of connected components (called "cuts" ) of the support of the densities of the spectra of the two matrices; in this 
terminology the case of genus spectral curves corresponds to the one cut assumption, whereas the higher genus 
case to multicut solutions. 

In the present paper we complement our previous paper [S] by computing explicitly the second and third order 
derivatives of the free energy. As we have anticipated, however, these result are derived independently of [S], only 
the notation being borrowed thence. We also point out that the resulting formulas are very close to similar residue 
formulas in |31| which, however, deals with a different hierarchy and moduli space and whose proofs are not given 
explicitly therein. 

The paper is organized as follows: in Section |21 we set up the problem and notation, recalling the formula 5 
for the free energy J- g over a spectral curve of genus g. We also link our present results to the current relevant 
literature. In Section 12.11 we give a to-the-bone review of the Bergman kernel on algebraic curves and of the 
properties that will be used in the sequel. In Section l2~2l we compute all second order observables, which include 
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the derivatives with respect to the filling fraction and the "temperature" t (which can be interpreted as number 
operator as well depending on the points of view) . Finally, in Section [5] we compute all third order observables 
by deriving a formula for the variations of the conformal structure of the spectral curve under variation of the 
parameters of the problem. 



2 Setting and notations 

As in our previous paper [3] we will work with the following data: a (smooth) curve S g of genus g with two marked 
points ooq, oop and two functions P and Q which have the following pole structure; 

1. The function Q has a simple pole at ooq and a pole of degree di at oop. 

2. The function P has a simple pole at oop and a pole of degree d\ at ooq. 

From these data it would follow that P, Q satisfy a polynomial relation and hence that the curve is a plane algebraic 
(singular) curve, but we will not need it for our computations 3 . Moreover we will fix a symplectic basis of cycles 
{dj, &j}_y=i,... s in the homology of the curve. By their definition we have 



di+l 



P = 



^ u K Q K - x - - - KUkQ K1 =■ V{{Q) - 1 + <D(Q- 2 ) , near oo c 



K=l ^ K=-l ^ 

d-2 + 1 — oo 



« = V ^ PK ~ 1 - p - E ^z^^ 1 =: V&P) -p+ 0{P- 2 ) , near oo P . (2-1) 

j=i j=-i 

The fact that the coefficient of the power Q~ x or P^ 1 is the same follows immediately from computing the sum of 
the residues of PdQ (or QdP) . The polynomials V\ and Vi defined by the above formula l|2-l|) are the potentials 
of the matrix model whose free energy we are considering 4 and the functions Q and P represent the semiclassical 
(commutative as per eq. 11-711 limit of the multiplication operators for the orthogonal polynomials. The coefficients 
UK, vj, t are read off eqs. (|2-lll as follows 

u K = - res PQ~ K dQ , vj = - res QP~ J dP , t = res PdQ = res QdP . (2-3) 

OOq OOp OOq OOp 

Note that the requirement that the curve possesses two meromorphic functions with this pole structure imposes 
strong constraints on the moduli of the curve itself. In fact a Riemann-Roch argument shows that the moduli space 
of these data is of dimension d\ + + 3 + g; to the above (di + 1) + (c?2 + 1) + 1 parameters we add the following 
period integrals referred to as filling fractions 



1 

2iTT 



PdQ, i = l,...g. (2-4) 



Here we have introduced a symplectic basis {a^, &j}i=i... g in the homology of the curve E ff and the choice of the 
a-cycles over the 6-cycles is purely conventional. 

3 We could also generalize this setting to the case di = di = oo which would simply mean that the two functions have one simple 
pole and one essential singularity; this generalization is quite straightforward and would need only some care about convergence. 
4 Formula 12-11 defines only their derivative; what we mean by V\(x) and V2(y) is explicitly 

v *v> ■= d £ u -if xK > v ^ ■= 15 ~T yJ ■ ( 2 - 2 ) 

K=l J=l 

The constant term in the two potential would have trivial consequences as it amounts to a rescaling of the partition function of the 
matrix model by e 11 o+"o j which of course has absolutely no consequences on all the observables of the model. 
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The free energy T g of the model (where the subscript refers to the genus of the spectral curve, not the genus of 
the h expansion) is then defined by the equations 

du K T a = 4 res PQ K dQ = U K , d v ,T g = \ res QP J dP = Vj , (2-5) 

K ooq 3 J OOp 



d €l T g = & PdQ =: T,. (2-6) 

Jbj 

These equations and in particular the derivatives with respect to the filling fractions appear in this precise context 
in |26j : nevertheless, after identifying the free energy with some tau function for a dispersionless Whitham hierarchy, 
similar formulas appear in |31j . 

The above equations mean that the coefficients of singular parts of P qua function of Q (or vice-versa) are 
the independent coordinates and the coefficients of the regular part (starting from Q -2 ) are the corresponding 
derivatives of the free energy. Moreover the a-periods of PdQ are the independent variables and the 6-periods the 
corresponding derivatives of T g . In fact one should add the extra constraints I\ = to this functional in order to 
assure that it comes from a saddle-point integration of the two-matrix model. However, since we are interested 
here mostly with the formal N — > oo, h — > 0, hN = t limit (as explained in (T^j) we do not impose the extra 
constraint and treat the filling fractions as independent coordinates. 

It was shown in (but see also ^5] for the extension to the chain of matrices) that J- g is given by the following 
equivalent formulas 



2T g 



I 9 

res P$idQ + res Q$ 2 dP + - res P 2 QdQ + i/x + V eT\ = (2-7) 

OOQ OOP 2 OOQ ^ ' 

i—1 

1 9 

K J IO ° Q (=1 

2 I\ 2 j ^ ^ 

K J i=l 

- » PdQ = -Vi(Q(0) + tln(Q(0) + / P&Q = u kQ~ K = OiQ-^-lO) 

J X i K=l 

-jQdP=-y 2 (P(C))+tln(P(C))+ y^QdP=^V: 7 p- J = 0(p- 1 ) (2-11) 









I In 


2l7T , 


/OOQ 








l In 


2iTT , 


/ OOp 



Q(0 
p 



j=i 
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p(0J 

li := res [Vi(Q) - t ln(Q/A)]d5 - res \V 2 (P) ~ t ln(PA)l dS - res PQdS" + V a I dS , (2-12) 

OOQ OOp L OOQ ^— ' J^. 

where dS* is the normalized differential of the third kind with poles at oopg and residues res dS* = — I = — res dS 

OOQ OOp 

and the function A is the following function (defined up to a multiplicative constant which drops out of the above 
formula) on the universal covering of the curve with a simple zero at ooq and a simple pole at cop 



A := exp dS) . (2-13) 



Remark 2.1 The "chemical potential" fi in the context of the dispersionless Toda hierarchy is the long-wave limit of the 
Toda field |4U| (denoted by cj> there) and satisfies, among others, the Toda field equation which is written, in our notations, 

' + dt exp(d t pO = . (2-14) 



duidv 



Remark 2.2 In formulas 12-71 . . . 12-121 the open paths of integration are supposed to not intersect the basis of cycles, i.e., 
to remain in the fundamental cell of the universal covering of the curve and the residues involving the multivalued function 
A are taken to be on the same cell of the universal cover. 
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We now make the necessary connections with the relevant literature |291 131)1 1321 1451 1471 1411 1371 • 

The equations that define here the free energy in the genus g = case are precisely the same that define the 
T-function of the dToda hierarchy where one imposes the (compatible) constraint of the string equation to the Lax 
functions. In the relevant literature |37l 1381 l3~9l 142| the functions P, Q are the Lax operators denoted by L, L^ 1 or 
C, C^ 1 and the normalization is slightly different. In higher genus our free energy is related to the tau function for 
solutions obtained via a Whitham averaging method on some invariant submanifold of the hierarchy. 

The link with the important works |32l 1291 1431 1471 141| is as follows; if the two potentials are complex conjugate 
Vi = V2 ='■ V then (in genus 0) the functions Q and P are conjugate-Schwarz-inverted in the sense explained 
presently; pick the uniformizing parameter A of the rational curve So to have a zero at oop and a pole at cx)Q(and 
a suitable normalization) 5 . Then one has 

Q(A) = PIT 1 ) (2-15) 

The function Q(A) is then the uniformizing map of a Jordan curve T in the Q-plane (at least for suitable ranges of 
the parameters) which is defined by either of the following relations 



Q(A) = P(l/A) or |A| = 1 . (2-16) 

In the setting of 0131213 E2] the function Q is denoted by z (and A by vS) so that then P is the Schwartz function 
of the curve L, defined by 

z = s{z) , z e r . (2-17) 

The coefficients of the potential V(x) = ^ are the so-called "exterior harmonic moments" of the region T> 
enclosed by the curve T 

t K = — I zz~ K dz (2-18) 

2l7T J T 

t = t = -J— [[ dz A dz = -J- I zdz = Area ^ , ( 2 -19) 

By writing ~z = S(z(w)) these integrals become residues in the w-plane. For conformal maps (i.e. Jordan curves) 
the r-function has been defined in |41j and given an appealing interpretation as (exponential of the Legendre 
transform of) the electrostatic potential of a uniform 2-dimensional distribution of charge in D [2H] 



ln(rr) 







i 


1 






z 


~ z 1 



d 2 zdV . (2-20) 



It can be rewritten as a (formal) series in the exterior and interior moments as (note that we are changing the 
normalization used in |291 1321 1451 14"7| to match more closely ours) 

21n(r r ) = -^// d 2 z\z\ 2 + t w + V (t K w K + t K w K ) (2-21) 



47T././z. 



K>0 



where the interior moments arc defined by (the normalization here differs slightly from |45j ) 



dt k rr =w K = -L JJ v zK <^ z . ^>0 (2-22) 
5 t0 r r =w a = ±-JJ^\n \z\ 2 d 2 z . (2-23) 

The logarithmic moment wq corresponds to our "chemical potential" //; it is not at all obvious (at least we failed to 
prove it by direct integration) but it is necessarily so because both J-q and it satisfy the same differential equations 
w.r.t. tK (or, in our notation uk = vk) and the same scaling constraint 

4^ = -t 2 + htdt + - k ) u k9u k + 5^(2 - J)vjd v Ar . (2-24) 

^ K J ' 



5 This A is exactly the same A appearing in the higher genus formulas; in fact —quite obviously— A = exp J 4^, which is translated 

A 



in the higher genus setting simply by replacing 4^ with the third-kind differential dS 



G 



Remark 2.3 Note that formula 12-71 1 is quite effective in that it allows explicit computations. For instance in the genus 
zero case one could do the exercise (e.g. with Maple) of explicitly presenting the free energy for arbitrary degrees of the 
potentials. In the approach of |29l 1321 1451 I47| the explicit computation in case of regions with finite number of nonzero 
external moments would be prevented mainly by the analytic computation of the logarithmic moment which, in our approach 
is nothing but fi — (Vi(Q) + V2{P)—PQ)o — tlri'y 2 , where the subscript o means the constant part in A and 7 is the "conformal 
radius" of the domain. This observation leads to the following 

Corollary 2.1 Let T> be a simply connected domain with finite number (say d + 1) of nonzer external harmonic moments 
and let 

d d 

z(w) = "/w + V 1 otjW~ j , S(w) = — + y aj-w 4 (2-25) 

3=0 3=0 

be the Riemann uniformizing map and Schwarz-function respectively. Then 

^JJ \n\z\ 2 d 2 z = (V'(z) + V 7 (S) ~ zS) a 



tln 7 2 , (2-26) 



where the subscript means the constant part of the Laurent polynomial in the bracket, t = ^Area(T>) and V(z) = 
Sk=i ^k~ zR '> w ith tx the K-th external harmonic moment. 

Recently J3U] the above tau-function Tp has been generalized to multiply connected domains. Also this case is 
contained in our case by imposing the same constraint on the potentials Vi = V2 = V. This amounts to saying that 
the curve admits an antiholomorphic involution ip given (coarsely speaking) by Q — > P and interchanging the poles 
oop and ooq. Such a curve would then become the Schottky double used in [j23J and the boundary of the multiply 
connected domain would be defined implicitly by Q(() = P(() (which defines -in certain cases- g + 1 cycles on the 
curve S g ) 6 . The reader familiar with the results in |30| will have no difficulty in finding the necessary dictionary 
to translate the two settings, provided he/she imposes the restriction on the potentials in our more general case. 



2.1 Bergman kernel 

The Bergman kernel is a classical object in complex geometry and can be represented in terms of prime forms and 
Theta functions. In fact we will not need any such sophistication because we are going to use only its fundamental 
properties (that uniquely determine it). The Bergman kernel f2(£, £') (where £, £' denote here and in the following 
abstract points on the curve) is a bi-diffcrential on S 9 x S s with the properties 



Symmetry: 
Normalization: 



n(c,c') = n(C / ,0 



C'ebj 



C') = 2iirujj(() = the holomorphic normalized Abelian differential 



(2-27) 
(2-28) 

(2-29) 



It is holomorphic everywhere on S 9 x S 9 \ A, and it has a double pole on the diagonal A := {£ = ('}: namely, if 
z(Q is any coordinate, we have 



(z(0-4C)) 2 



So(0+0(z(0-z(C')) 



dz(C)dz(C) , 



(2-30) 



where the very important quantity <So(C) is the "projective connection" (it transforms like the Schwarzian derivative 
under changes of coordinates). 

6 An algebraic curve of genus g with an antiholomorphic involution and g + 1 real components is called an Harnack's M-curve; in 
general a curve with the aforementioned properties can have k < g + 1 real components, see for instance !24l 

7 Our use of the term "Bergman kernel" is slightly unconventional, since more commonly the Bergman kernel is a reproducing kernel 
in the L 2 space of holomorphic one-forms. The kernel that we here name "Bergman" is sometimes referred to as the "fundamental 
symmetric bidifferential" . We borrow the (ab)use of the name "Bergman" from 1271 . 
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It follows also from the general theory that any normalized Abelian differential of the third kind with simple 
poles at two points z_ and z + with residues respectively ±1 is obtained from the Bergman kernel as 



ds z+ , z _(o= n(cc')- 



(2-31) 



C'=a- 



2.1.1 Prime form 



For the sake of completeness and comparison with the results of |3()j we recall here that the definition of the prime 
form £(C,C')- 

Definition 2.1 The prime form is the (—1/2,-1/2) bi- differential on E g x S g 



E((,C) = 



\ a 


(u(C)-u(C')) 


VP 


](0^-](C) 



fc=l 







A 





where LOk are the normalized Abelian holomorphic differentials, u is f/ie corresponding Abel map and 
half-integer odd characteristic (the prime form does not depend on which one). 

Then the relation with the Bergman kernel is the following 



(2-32) 
(2-33) 
is a 



tt(0C')=d c d c ,ln£(C,O= E 9 »A 

k,j=l 



a 




A 





^(CH(C') 



u(C)-u(C) 



(2-34) 



2.2 Second order observables 

The second order observables of the two-matrix model (in fact, the multi-matrix model) have already been inves- 
tigated in the literature ^3 [3U1 an d their relation with the Bergman kernel extensively documented. Here we 
just bring a different and possibly more rigorous derivation of those identities. The main reason why the second 
order observables appear already in the literature is the expectation (and -in some cases, mostly in the one-matrix 
model setting- proof) of their "universality", that is their "independence" on the fine details of the potentials. This 
paragraph will support once more this point of view; indeed we will see that these generalized "specific heats" do 
not really depend on the potentials but only on the spectral curve of the model and can be described by geometrical 
objects directly linked only to the curve itself. That is to say that these second order observables will be the same 
for any pairs of potentials for which the conformal structure of the spectral curve is isomorphic. 

Of course this is no proof of universality, since there is no (rigorous, to our knowledge) proof that the free energy 
as defined in this paper is really obtained from some scaling limit of the partition function of the matrix model. 
Such a proof exists in the one-matrix model and uses rather sophisticated tools (Riemann Hilbert problem) |15j . 
A first step in this direction has already been taken by collaborators and the author in [7] . 

On the other hand we will see in Section |21 that universality does not hold for third order observables, which 
will not depend purely on the conformal structure of the spectral curve but on the functions P and Q explicitly 



Let us start with d UK d Uj T 



2i-Kd UK d Uj T = — d u , 



PQ dQ 



K<t> {du,P) Q Q K dQ 



(2-35) 



Let us focus our attention on the differential {d UK P)c)dQ (or (d Vj P)QdQ), where the subscript indicates that the 
corresponding quantity is kept fixed under differentiation. It follows from the definition of the coordinates on the 
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moduli space and from eqs. (|2-l|l that 

{d UK P) Q dQ = (Q K - 1 + 0(Q- 2 )) dQ near oo Q (2-36) 
(d UK P) Q dQ = -(d UK Q) P dP = 0(p- 2 )dP near oop (2-37) 

I (d UK P) Q dQ = d UK e, = , (2-38) 

and 

(d Vj P) Q dQ = -{d Vj Q) P dP = (-P' 7 - 1 + 0{P- 2 )) dP near oo P (2-39) 

(d VJ Q) P dP = 0{Q- 2 )dQ near oo Q (2-40) 

<f (d VJ P) Q dQ = -d VJ e j = . (2-41) 
J a,j 

In these formulas we have used repeatedly (and we will use it many times) the so-called thermodynamic identity 
(or reciprocity) 

(dP) Q dQ = -(dQ)pdP , (2-42) 

where d denotes any derivative. This is immediately obtained by differentiating P as a composite function with 
the (local) inverse of Q. 

Note that the only singularities of these differentials are at the two marked points (see the discussion on the 
third-order observables for a proof that there are no poles at the branch-points). In other words {d UK P)QdQ and 
(d v ,P)qdQ are Abelian differentials of the second kind (i.e. with poles but no residues), normalized (i.e. the 
a-periods vanish) with a pole of degree K + 1 or J + 1 respectively only at ooq or oop respectively. It is immediate 
that this differential is uniquely determined and can be expressed in terms of the Bergman kernel as follows 

(d UK p) Q dQ(0 = — L- I Q K ((MC, , 

LITT J\ J ooq 

(d Vj P) Q dQ(() = -(d Vj Q) P dP = j P J (CMC, C) ■ (2-43) 

Indeed, when £ ~ ooq we have (using z = 1/Q(£) as local coordinate) 

1 f _r„,. 1 



/ z-'ntC, C) = -^— (fz- J 1 + 0(1) dzdz' = (-«- J - 1 -H9(l)) dz = (Q J -^(Q- 2 )) dQ (2-44) 
ZiirJ J Hit J J \_(z z ) 

It then follows that 

d Uj d UK F= 1 / / Q K (C)Q J (CMC,C) ■ (2-45) 

[ZITT) J\J J ooq J ooq 

By similar arguments one obtains also the formulas 

du K d VJ T Q K (d Vj P) Q dQ t ~^- im (2-46) 

Z17TA Joo Q 

y Q K {d Vj Q) P dP = (2-47) 

tt(C,C)Q K (C)^ 7 (C) (2-48) 



2iirK . 
-1 



{2iir) 2 KJ 



OOQ J OOp 



We now compute the second derivatives w.r.t. the filling fractions; to that extent we notice that (d <Lj P)QdQ is 
holomorphic everywhere because 

(8 PWD - / °(Q~ 2 ) d Q near °°Q fo 

(d ej F) Q dg - | = _p Q)pdP = 0(P -2 )dP near ^ . (2-^) 
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Thus it is regular at the marked points and has no other singularities (the poles of the variation at Q fixed of P 
cancel with the zeroes of dQ; see discussion in section |3J). Moreover it satisfies 

2in6 jk = 2md^ k = d t] I PdQ = I (d ej P) Q dQ . (2-51) 

J a k Ja k 

Therefore (d tj P)qdQ = 2iiuxij where {^j}j=i...g is the basis of normalized Abelian holomorphic differentials. As a 
consequence we obtain 

d t] d tk T = I (d £j P) Q dQ = 2wr I ojj = I I n = 2mr ik , (2-52) 
Jb k Jb k Jb k Jbj 

where r = [T%j]i t j=i... g is the period matrix for the holomorphic curve. 

The other mixed derivatives are easily computed along the same lines as above to be 

dui{de ^ ^ 2^K f ceb f n((X')Q K (C) (2-53) 

dA ^=2^Ki Q QKdS (2 " 54) 

dtd VJ T = / P J dS (2-55) 

d t d Cj T = f dS = Uj-(ooq) - Uj(oop) (2-56) 

Jbj 

9 t 2 ^-ln( 7 7), (2-57) 

ln( 7 ) := - res In ( ^ ) dS , ln( 7 ) := res In (FA) dS . (2-58) 

ooq \ A J oop 

(Note that the product 77 does not depend on the arbitrary multiplicative constant entering the definition of the 
multivalued function A). The first two relations are obvious recalling that |S] 

dS = -(d t P) Q dQ = (d t Q)pdP . (2-59) 

The third relation is a well known result for the third-kind Abelian differential. 

We shall prove the last relation: we start from the first derivative of T (/i, also called chemical potential) and its 
expression given in our previous paper |S] 

2i-n\i = 2md i ,T=l (Vi(Q) - t In (Q/A)) dS - I (V 2 (P) — t\n(PX)) dS — I PQdS + 2m ^ ej IdS = (2-60) 

•/ oon J oop J oon „• 1 J 



Q j=l ■> bj 



= I In (Q\ PdQ + I In (PA) QdP - <f PQdS + 2m^t il <f dS 
We can therefore compute 

: — 2i7T In 7 



(2-61) 



2ind t fi = -<b \{d t \) Q PdQ+ <f> \n(Q\{d t P) Q dQ+ (2-62) 



ooq n Jcoq \ 



—2iir In 7 



=dS 

j(d t X) P QdP+ I ln(PA)'(^g)7d?+ (2-63) 

OO p ™ J QQp 
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<9t*> 



-dQ 



-I P (d t Q)x^ ~ I Q {d t P) x ^ +2^^ <f (d t dS) P = 

Jooq " J ooq " j =1 Jbj 

= 2in 111(77) j 

where we have used Riemann bilinear identity to move the residue of the following term on line fHty 



(2-64) 
(2-65) 



— — 2?7T£j 



/ \ (d t \) P QdP = - I j(dtX) P QdP + ^2<f QdP 

J OOp ^ J OOp ™ j — \^ a 3 



<f (d t dS)p - <f QdP<f (d t dS)i 

Jbj Jb; J a,j 



=d t { dS=0 



(2-66) 



This concludes our analysis of the second order observable of the model. We repeat here for the sake of completeness 
that we are aware of a prior derivation in |19| and |46) for the case of conformal maps and genus (one-cut case 

Hi)- 

2.3 Extension to formal power series for the potentials 

The present derivation for the second derivatives does not rely on formal manipulations involving higher uk or vj. 
However it gives a deeper insight to use such a formalism. To that end we will now think of the two potentials 
V\ and V 2 as infinite power series, without concern about their convergence. This is often done in the physical 
literature 02] and yields the correct results in a faster and possibly more elegant way. To this end we define 
the puncture operators 



E t*- 1 

K=l 



Kd UK 



5V 2 {p) 



,7=1 



Jd Uj 



(2-67) 



Using this formalism we realize that our free energy T g is nothing but the generating function of the Bergman 
kernel on the selected curve 



5 2 T 



SV 1 (q)SV 1 {q') 



1 

"4^2 



1 



fi(C(<?U(g')) 



(q - Q(0W - Q(O) (q-q') 2 dQ(C(<z))dQ(C(g')) 



(2-68) 



where the integrals are to be read in a formal sense of inverse power series in q, q' , coefficient by coefficient; indeed 
each coefficient is precisely given by formula ()2-45JI . In eq. (|2-68|l the notation £(g) means the point on the curve 
that projects on the physical sheet of the projection Q. The physical sheet is the sheet that has no cuts extending 
to infinity. On the physical sheet the formula is not formal; on the other sheets it should be properly understood 
as analytic continuation. 

Similarly we have the other formulas where the P-projection is involved 



5 2 T 



5V 2 (p)5V 2 ( P ' 
5 2 T 



jrdpdp' = --^2 dp dp' I I 

/ J OOp J OOP 



1 



SV 1 {q)5V 2 (p') 



-dqdp' = -^dqdp'f £ 



= n(C(p),C(pO) 

(p-p(C))(p'-p(C0) {p-p'Y dP(c(p))dP((V)' 
n(CC') n(c( q ),ap')) 



3-69) 



OOQ J OOp 



(q - Q(0)(p' - P(O) dQ(C(g))dP(C(p')) 



(2-70) 



To be more precise, these formulas are exact (modulo the issue of convergence of the potentials) only in a neigh- 
borhood of the respective punctures where Q and P provide univalued coordinates. 

It is interesting to recast the problem upside-down: let it be given a smooth curve of genus g S 3 with Bergman 
kernel f2 (which implies a fixation of symplectic basis in the homology of the curve) and two marked points ooq^p 
such that: 
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there exist two functions Q, P with simple poles at the respective marked points; 

the only other singularity for the functions above are essential singularities at the other (respective) marked point. 

In this setting we can define "coordinates" uk , vj by means of formulas (|2-3[1 and then the expansion of the 
Bergman kernel around the punctures in the coordinates provided by the functions P, Q are the second derivatives 
of a "tau" function. 



2.3.1 Relation with conformal maps 

The reader trying to bridge the present setting with should be warned that the puncture operator here is not 
the "infinitesimal deformation" operator ibidem. Their operator would be written in our notation 

oo e oo 

E ■ ^7773 = E «TX, , (2-71) 



SV{(x) ^ UK > 5V'(y) f[ 



with the relation (e.g. for the x — x two-puncture case) 

d d S 5 r S 6 ^ , nfm . 

d Xl dx 2 5V{( Xl ) SV{(X 2 ) 9 SV^xt) 5V l {x 2 ) 9 [ ' 

Therefore their double-deformation is the antiderivative of our double puncture. This is why in |3(J| they have the 
logarithm of the prime form 8 whereas we have the Bergman kernel. Since in this paper we will not investigate 
Hirota equations (i.e. Fay's identities) which involve the prime form more naturally than the Bergman kernel, we 
will not pursue this observation any further. 



3 Third derivatives: Beltrami differentials 

In this section we compute the third derivatives in terms of the canonical structures and data of our moduli space. 
We must note here that in the case of the tau-function for conformal maps of connected domains (i.e. genus 
case) constructed by Marshakov, Krichever, Wiegmann, Zabrodin et. al. these formulas are known [22] (see 
appendix [5J|. In order to obtain their "residue formulas" from the formulas we are going to write down, we should 
restrict ourselves to the genus zero case and impose that P(A) = Q(l/A), A being the uniformizing parameter of 
the rational curve with A(ooq) = oo and A(oop) = (see appendix lB|l . However our formulas cover a more general 
case that would includes the situation of conformal maps of multiply-connected domains |30| . were we to impose 
some restriction of reality on P and Q (as explained in section EJ) 

When computing a variation in the moduli of our problem {uk, vj, t, ej}, we introduce a deformation of the 
conformal structure of the curve. We will see that this deformation is equivalent to introducing simple poles at 
the branch-points of one or the other of the two projections Q : S g — > CP 1 and P : S g — > CP 1 , depending 

on which among Q and P we keep fixed under the variation. Let us denote by q^, /i = 1, , d% + 1 + 2g and 

p v , v = 1, , d% + 1 + 2g such points and then by Q^, P v the corresponding critical values; here we are implicitly 

making an assumption of genericity of the potentials, in that we will consider these branch-points as simple. Were 
we to consider more degenerate cases (critical potentials) we would have to modify some of the formulas below 
(mainly the Beltrami differentials). Note that critical potentials (i.e. potentials which are fine-tuned so as to have 
degenerate singularities of the maps P and Q) are very important for applications to conformal field theories in 
that they provide instances of minimal conformal models. 

Coming back to our generic situation, let us consider the case of the Q-projection, the other case being totally 
similar. The projection Q : S ff — ► CP 1 determines the conformal structure of the curve and hence it suffices to 
consider the variation of the structure arising from the deformation of this projection. Consider first the branch- 
points of the Q-projection and introduce the local parameter in a neighborhood of the critical point 

^ = VO~Q^ ■ (3-1) 

8 Here we are speaking rather loosely; we should also fix a coordinate and consider the prime form as a bi-function by trivializing 
the spinor bundle. 
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As we know, the critical values of both projections Q and P are not independent coordinates in our moduli space; 
they vary when performing a variation of our moduli. Let us call d any infinitesimal variation of our moduli. 
Under such a variation Q M also undergoes some variation dQ^. Suppose now we have a germ of function F in 
a neighborhood of (which may also depend explicitly on the moduli) and we make the variation d at Q- value 
fixed. Then 

(dF(C,)h = + d d F(C,) , (3-2) 

where doF denotes the variation of F coming from the explicit dependence on the moduli. We see immediately 
from (|3-2[1 that, generically, the variation has a simple pole at each branch-point of the Q-projection. Ditto for the 
branch-points of P when computing a variation at P-value fixed. 

A variation of the moduli induces a variation in the conformal structure of the underlying real surface E s ; these 
infinitesimal variations are described by the so-called "Beltrami differentials" . They consist of differentials of the 
form 

/j,(z,z)dz/dz e £°°(E 9 ) , (3-3) 
and enter into Rauch's formula I35| for the variation of the Bergman kernel 

W,C') = / / s KC'MCCMCC) ■ (3-4) 

It is shown in many places |221 1351 1771 128| that varying the image of a branch-point b of a covering ip : S ff — > CP 1 
corresponds to the following Beltrami differential (in terms of a local coordinate z centered at b) 

VSc = , (3-5) 

where l e is the characteristic function of the e disc and e is a parameter small enough so that no other branch-points 
fall within the disc. This is called "Schiffer" variation and corresponds to varying the critical value of a simple 
branch-point; for higher order branch-point see |27j . Plugging this into l)3-4|l and using Green's theorem followed 
by Cauchy's residues theorem one obtains simply a residue 



n(CX"MC X") 
T d P (C") 



<W"(C, C) = res • (3-6) 



Our setting is slightly different, in that we have two projections P and Q and the respective branching values 
are not independent moduli but vary together with the {uk, Vj, e^, t} moduli of our problem. Therefore all we need 
in order to find the Beltrami differentials corresponding to these variations is to be able to find the coefficients dQ^ 
and dP v that appear in eq. (|3-2|l . To this end let us now consider the case of a differential of the form FdQ. It 
follows from the above discussion (and we have already used this fact implicitly to compute the first variations of 
PdQ) that its variation at Q fixed will be holomorphic at the branch-points because the differential dQ has simple 
zeroes there which cancel the poles of the variation of F. 

Suppose now that for some reason we have a way of identifying independently what differential is to := (dF)gdQ 
Then we could compute the desired coefficients 9Q M as in the following formula 

u dF 
3^Wa») = -dQf. — iq^) that is... (3-7) 

dQ» = -^p{%) (3-8) 
We now apply the above to the case F = P and the differential PdQ; in this case we have already discovered that 

1. the derivative of PdQ w.r.t. ej is the holomorphic normalized Abelian differential (d tj P)qdQ — u)j\ 

2. the derivatives of PdQ w.r.t. uk or vj are given by eq. (|2-43|l : 
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3. the derivative of PdQ w.r.t. t is the normalized Abelian integral (<9 t P)gd<2 = dS of the third kind with poles 
at the marked points. 

Therefore the parameter dQ^ and dP v are given by the following formulas 

n(C,C) 



dP 



1 



1 Jceb 



q . 2m dP(C) 



2i 7 rii' dP(C) 



9 t Qu = - 



2inJ dP(C) 

/■OOp 

/ "(00 

J OOQ 



dP(C) 



(3-9) 



(3-10) 



(3-H) 



(3-12) 



Repeating all of the above but interchanging the role of P and Q we find (note the opposite signs, coming from 
the thermodynamic identity) 



BP — 
de > u ~ dQ 



1 



"(CO 



r) P — 



1 



2in dQ(C) 

Q(() K n(cc) 



2mK dQ(C) 

P(C) J 0(CC) 



a p - - 



d t P, 



2i7rJ dQ(C) 

""(00 



dQ(C) 



C=p, 



In general we can write the formula (d denoting an arbitrary variation) 

(0P) Q dQ 



3Q M = - 
dP y = - 



dP 

(dQ)pdP 



in 



dQ 



= (8P)flW 



(3-13) 
(3-14) 
(3-15) 
(3-16) 

(3-17) 
(3-18) 



We are now in a position of expressing the Beltrami differentials corresponding to a variation in our moduli as in 

(3-19) 



d 2 + l+2g di + l+2g 

m= E 9Q^ } = E 9p v ^, 



v=l 



where the notation jj,^ denotes the Schiffer variation centered at the point (. The equivalence sign = means 
that the two expressions determine the same variation of conformal structure. Indeed we should remind that any 
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two Beltrami differentials /i and fi determine the same variation of the conformal structure if for any quadratic 
holomorphic differential <E> we have 

" /i$ = JJ . (3-20) 
Indeed, for $ S H°(K?) a holomorphic quadratic differential we have 

d 2 + l+2g d 2 + l+g d 2 + l+g , „ , 

til JJz. c-i.de} ^ c=«„ dQdP 

- § S 2 ^^ E E <-) 

thus proving the equivalence of the two Beltrami differentials. The equality marked (■*•) follows from the fact that 
the differential ^jqgjr - has poles only at the branch-points of P and Q, and the sum of all residues must vanish 9 . 

From Rauch variational formula H3-4|) (see also |27| ) follows then that the variation of the Bergman kernel under 
such a Beltrami differential is 

(«kc,o) - j e% <».^p- - t.» (ap)g z°^r <c, ' a <- 3 » 

^ C=<z„ dQ(C) ^ C=<r„ dP(C)dQ(C) 

^ C=p, dP(C) t^i C=p» dP(C)dQ(C) 

Using the expressions above for the coefficients 9Q M and 9Pj/ we obtain then 



, d 2 + l+2g f fi(C,CMC,C"MC>C') 

(a e< n(C,O)0 = y res = = (3-25) 

3 9 2 ^ ^ dP(C)dQ(C) 

d 2 + l+2g 

(a „ n(CO)0 ._ g dP(?W(?) <-) 

i 2 +i+2 9 { p(C') J n(cx)n(lc"MlC) 



(d Vj fl(CX'))o = 7^— T Y res J ("~°° p ^ = (3-27) 

' w ^ f= ,„ dP(c)dQ(o 

fi(C,0"(C,C"MC,O 

t^i C=«m dP(C)dQ(C) 

The variations at P fixed would be given by the same formulas as above but replacing the sum over the critical 
point of Q by the sum over the critical points of P and changing the overall sign. 

The reason why we distinguish between the variations at Q or P fixed will be clearer in a moment. 
Indeed let us compute now the third derivatives of T g . If for example we want to compute d UJ d UL d VK J- g the 
simplest way is to leave the variation of vk last so that we will be varying at Q fixed; 

3uj3ul&vkJ~ g — ^vk^uj^ul^~ g (3-29) 

= ^~d VK I 9l ( d UL P) Q dQ = —±—-d v J I Q(0 J Q(C) L n((X') = (3-30) 



9 The differentials (dP)QdQ have poles at most of degree d,2 + 2 at ooq or di + 2 at oop but those are canceled by the poles of 
dPdQ at those points, of degree <fe + 3 and d\ + 3 respectively. 
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' ' ' Q(0 J Q(C) L (d VK n(( 1 C)h= (3-3i) 



-^T - I / / (3-32) 
8in s JLK ^ c=?, /oo Q Joo Q Joop dP(C)dQ(C) 

d2 ^ 29 (9^P)odg(^ J P) Q dQ(9 I1 ,P) Q dQ 



Performing the variation at P fixed, since P is not a local coordinate near the marked point ooq, would introduce 
unnecessary complication to the formula (basically adding an extra residue to it). Thus we need to pay the necessary 
care in choosing which among P and Q needs to be kept fixed in order to have the simplest formulas. 

With this caveat in mind can thus now compute the third derivatives of the free energy. The resulting formulas 
should be compared with the residue formulas in EH| (which would correspond to genus g = here) The 
formulas are simplified if we introduce the quadri-differential 

0(3,1),, > f . n _ «MM f ,-ui 
n (Ci,C 2 ,C,0- d p(C)dQ(o ' (3 " 34) 



and the notations 



'ooq 



2iirK 
1 



We also recall that 



Vj(C) := (d Vj P) Q dQ(0 = <f> P(C') J tt(C,C) • (3-36) 

2nrJ J OCp 

= (d e] P) Q dQ , dS = -(d t P) Q dQ , (3-37) 

are the normalized Abelian differentials of the first and third kind respectively. With these definitions and reminders 
we have 

- T fl y -^. (3-40) 



1 ti2 + l+2g .oop 

0*0^* = -7^2 E , res ? f / ^ (3 ' 1) (Ci,C 2 ,C 3 ;C) = (3-41) 



3=°°Q 



d2 ™ ^(C)^-(C)d5(C) _ dl ^ 29 ^(C)^-(C)dg(C) 
^c=t dP(C)dQ(C) " ^c=p„ dP(c)d (C) ( '" 42) 



d2 + l+ 2g /.oop poop 

de^Tg = -— E res i / / n( 3 ' 1 )(Ci,Ca,Ca;0= (3-43) 



^ V^res ^» - ^V^res (3 44) 

2- c=i dP( C )dQ(C) - ^ c r 4l dP(C)dQ(C) {6 44) 
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The third derivative w.r.t. t has to be treated with some care, although the result is a formula of the same sort. 
For the purpose of this computation it is best to choose as basepoint defining A the marked point ooq, regularizing 
appropriately the integral 

lnA=/dS:= lim [ dS - ln(Q(0) = lim f dS - ±- ln(P«)) + -J- ln("«b+i) ■ 

J ooq C^°°QJ( C^oOQjf «1 «1 

This is just a convenient normalization of the function A with the property that 



(9 t ln(A)) P = f (d t dS)p , (9 t ln(A)) Q = f (d t dS) Q . 

J ooo J ooo 



'OOQ -'OOQ 

Using these properties we find 



d z t T g = d t ln( 7 7) = d t ( res ln(Q/A) — - res ln(PA) — 



res ln(PA) — 

OOQ A OOp A 

{d t Q)xd\ {d t P) x d\ 
res — n\ rcs p\ — = 

ooq Q\ oop p"A 

(d t X) Q dQ (d t X) P dP 

ooq Q\ oop f\ 

(9 t ln(A)) Q dQ (a t ln(A)) P dP 

- res + rcs = 

ooq y oop F 

- res — - h rcs — - 

ooq Q OOp f 



=0 



/■ooq roOQ roop roo P 

/ / (d t Q) Q ~ / / (9 t O) P - 

J OOQ J OOp J OOQ J OOQ 

pOOp r OO p fOOp 

E , res / / / n (3,1) (Ci,C2,C 3 ;C) = 

^=1 ^ Qfi •> Cl=OOQ J C2 = °°Q J C3=°°Q 

d g( C) 3 *^+2 fl d5(c)3 

2^ f=n.. dP(t\dO(t\ ~ 2s t V =f, 



The other derivatives are 



^ c=*, dP(C)dQ(C) ^ c= P , dP(C)dQ(C) 



d d d t - - V res ^i^h^w. = v 



rcs 



^ c=g „ dP(c)dg(c) ^ c=p« dP(C)dg(c) 
hl+2 " ^(C)^(C)V,(C) _ dl ^ 2ff _ «Vj(0 



9 9 3 T - - V res - V rc 



rcs 



^ c=*, dP(c)dQ(C) ^ c=p, dP(c)dg(c) 



9 a a ^ - - V res ^AW"iU^JV^ _ \ - 



rcs 



^ c= 9 „ dP(C)dg(c) ^ c= P , dP(c)dg(C) 



' dS(C)zMC)Vj(C) d1 ^ 29 d5(C)^(C)Vj(C) 



rcs 



^ c=* dP(c)dg(c) ^ c= P , dP(c)dg(c) 

In all of the above formulas we can move the sum over the residues at the critical points of cither maps P and 
g because the differential has poles only at those points and not at the marked points. Indeed the denominator 
has a pole of order d\ + 3 and c? 2 + 3 at ooq and oo p respectively, and the pole of the numerator never exceeds 
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these values. In the remaining derivatives, instead, the simplest form is obtained by summing over only the critical 
points of one of the two maps, case by case: 



^ T dS(Q*U K (Q 
1 UK 9 ' ^C=l dP(C)dQ(C) 

(J, — ± 

dtdvAs - h ™l dP(OdQ(o 



8 t 8 M = T 9 res 



^ C= 9 , dP(C)dQ(C) 

d d T - V^res "«(CWWO 

" UK UJ 9 ~ ^C=l dP( C )dQ(C) 

„ „ „ - _ d2 v^ 29 WW(C) 

Ou^d^J-g KS dP(C)dQ(C) 

dv^dujJ-g ^ K» d p (C)dQ(C) 



ft* a ^ - -y'^res ^(0^(0^(0 
<kd VK d VJ J- g 2^ ™l dP(C)dQ(C) 

W**"*' ~ ljc~ dP(C)dQ( C ) 

a a a ^ - dl V 29 res 
" l , " " J 9 ~ ^ f 4l dP(c)dg(0 

Wl(C)Vk(C)v./(C) 



d UL d VK d Vj T g = II 



res 



^ c= P „ dP(c)dg(0 



3-59) 
3-60) 

3-61) 
3-62) 
3-63) 
3-64) 

3-65) 
3-66) 
3-67) 
3-68) 



Although we have not written it each time, it is clear that any of these derivatives is obtained by means of suitable 
integrals of the kernel f^ 3,1 ^*, •, •; Q over 

1. a 6-cycle corresponding to each variation w.r.t. e's (times 7^7); 

2. a circle around ooq against Q K /(2in K) for each variation w.r.t. m^'s; 

3. a circle around oop against — P J /{2i-K J) for each variation w.r.t. uj's; 

4. a path from oop to ooq for each variation w.r.t. t; 

followed by minus the sum over the residues at the points ( — if there are two or more derivatives w.r.t the ur-'s 

or two derivatives w.r.t. t and one w.r.t. uk, 

or 

the sum of residues at the points C = lV s if there are two or more derivatives w.r.t. the vj's or two derivatives 
w.r.t. t and one w.r.t. vj. In all other cases the choice of residues is immaterial. 

Using the above results somewhat liberally in the case of potentials given by power series we obtain 
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Corollary 3.1 [Three-loops correlators] Denoting by the inverse function of Q(Q — x and £(y) the inverse 
function of P(£) = y on the respective physical sheets we can write the three-loop correlators as follows 



Tr- 



1 



Tr- 



Mi 
1 



-Tr- 



1 



Mi 



-Tr- 



Mi 
1 



-Tr- 



Mi 



conn 



M 



-Tr- 



Ms 



V res ^ 3 '%(x),C(z'),C(V');C) 

dg(c(x))dQ(c(x'))dg(c(x")) 

^ res fi^Cto.C^UMO 

dQ(c(x))dg(c(^))dP(C(^")) ' 



(3-69) 



(3-70) 



and similarly for the other two three-loop correlators but summing over p v 's and with an overall minus sign. 

"Proof". The quotes are because a complete proof should address also the issue of convergence of the series 
involved. Indeed we can obtain the above correlators by applying the relevant puncture operators, which can 
be done only by considering infinite series for the potentials. At the end of this process we should set all the 
coefficients but a finite number to zero so that the final formula involves only a finite sum of residues. Since we are 
not addressing the summability of the series involved in this procedure we cannot claim rigor for this proof. With 
this in mind we nevertheless proceed on a formal level. For instance the first correlator would be 



E 

K,K',K"=1 



x -K-i x f-K'-x x n-K- 



- x KK'K"d UK d UK ,d UK „T g 



(2inf 



E 



K,K> K"=\ 



y res 



Q(0 K Q(C) K Q(C') K 0(3 , 1)r > r „.,, M 

rK r x l\K'+l / X H\K"+1 I^jS, >W — W 



(3-71) 
(3-72) 



We can interchange the sums with the integrals (and sum the resulting geometric series) if the loops around ooq 
project -e.g.- to circles in the Q-plane which leave x, x' , x" in the unbounded region, that is, if \Q(Q/x\, etc. remain 
bounded and less than one. In this case we have then 



-1 



E 



i 



i 



i 



C=?m x - Q(C) x> - Q(C) x" - Q(C") 



(3-73) 



which provides the result by residue evaluation (the only residues are at the three points ((x), C( x '), C( x "))- The 
other correlators are computed exactly in the same way. "Q. E. D." 



4 Conclusion and outlook 

The main result of the paper is the formula for the triple derivatives of the free energy for arbitrary genus spectral 
curves. En route we have computed the variation of the conformal structure of the curve, which would in principle 
allow to compute derivatives of any order. Quite clearly the task becomes an exercise in complication (see similar 
problems in |46j ) . We also remark that in many practical uses the third derivatives suffices to determine many 
relevant properties. In particular our formulas should be sufficient to address completely the issue of associativity 
in both sets of variables {uk} and {vj}. It should be remarked that the multi-loop correlators of the one-matrix 
model (which is a sub-case of the two-matrix model when one potential is Gaussian) were presented by means of 
a recursive (and hence not totally explicit) procedure in 
We have not considered the mixed correlator 

M(x,y) := (tt \ \ ) . (4-1) 
\ x-Miy-M 2 / 

Such correlator is not universal and can be obtained via loop equations in the large TV- limit [361 121)1 1171 1191 . 
The mixed correlator l|4-l|l has also already been computed exactly in the finite N = dim(M^) regime in |10j . 

Acknowledgments. I would like to thank Bertrand Eynard, John Harnad, Jacques Hurtubise, Alexey 
Kokotov, Dmitri Korotkin, Alexander Bobenko for stimulating and critical discussions. 
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A Example: the Gaussian— Gaussian model 



If both potentials are Gaussian (d\ = d 2 = 1) then the spectral curve is a fortiori rational. There is nothing 
new in what we will write in this appendix since the tau-function has already been computed (see |29|. with the 
warning that we are using different normalizations for the parameters) corresponding to the elliptic case (it suffices 
to consider the holomorphic and antiholomorphic coordinates as independent in their formulas to obtain ours). We 
choose the normalization of the coordinate A so that 

Q(A) = 7A + a + a 1 X~ 1 (A-l) 
P(A) = 1 \- 1 + (3 + PiX , 7 eR+- (A-2) 

Using the formula for the free energy and computing the residues one obtains 

2T = tHn( T ) **» + f (W+^ + 2r, lUl ) 
V Vu 2 v 2 - 1 / 2 u 2 v 2 - 1 

In this case one can check by direct computations that our residue formulas yield correctly the third derivatives of 
To (a computer aided algebra system like Maple will help for the computations of residues). 



B Relation with conformal maps 

We want here to show that the "residue" formulas in |461 are contained in our formulas. We recall that the 
comparison is achieved by taking (see Section O z = Q, d\ = d 2 = d, V\ = V 2 = V and using 

d 

Q(A)= 7 A + ^a J A- i (B-l) 

i=0 
d 

P(X) = 1+J2^ (B-2) 

7 £ R+. (B-3) 

In this case the Bergman kernel is simply 

n(A,A') = ^ 5l (B-4) 

and the differentials of the second kind reduce to 

U K = (d UK P) Q dQ = 1 jes^niX, X')Q(X) K = (Q K (X)) + d\ = ±d(Q K (X)) + , (B-5) 

where the subscript ()+ means the nonnegative part in A. Then our formulas reproduce exactly those in |.'{2j when 
proper identifications are made (our normalization differs from theirs). 



C Example: the Gaussian— Cubic case 

Here we show how to compute the free energy and derivatives for the case V 2 (y) = -yj/ 2 + v\y and V\ cubic (one 
could do the exercise for higher degree potentials, the expressions become quite large but can be obtained explicitly 
e.g. with Maple). This case corresponds to the one-matrix model 10 and indeed the spectral curve is hyperelliptic. 

10 Quite clearly the corresponding matrix model would have to be defined with normal matrices with spectrum on certain contours, 
rather than Hermitian matrices 
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We will however restrict ourselves to the genus zero case where computations can be made completely explicit. 
Indeed we have 

Q(A)= 7 A + ao + y (C-l) 

2 

We see that here the computations of all derivatives are possible in explicit terms because Q has only two critical 
points; in all cases in which (when computing the third derivatives) we should sum over the residues at the critical 
points of P wc can replace the sum by summing over the critical points of Q plus the other residue at A = 0. Wc 
can eliminate the parameters at as follows (sometimes that makes formulas shorter) 

do = vi + v 2 (3o, oil = «27 i t = oiiPi ~ 7 2 ■ (C-3) 

In order to express the free energy in terms of the coordinates Uk we should invert an algebraic set of equations 
which can be done numerically at least. Nonetheless wc can compute the free energy and any derivative up to the 
third using our formulas. The quadri differential fi 3 ' 1 is simply 

fi( 3 ' 1 )(A 1 ,A 2 ,A 3 ;A) = (Ai _ A)2(A 2 dA A^(t A - d A)2p'(A)g'(A) (C " 4) 

Using e.g. Maple one can immediately compute all derivatives by residue evaluations. For instance the chemical 
potential is 

H = dtFo = 7 - "o/3o - Pioli + 2 1 In (7) + (C-5) 

27 37 

and we have (just two examples) 



dt ^ 7 5 - 4 /? 2 2 «2 3 7 3 + 7 3 «2 2 /3i 2 - 2 7 4 « 2 /?! (C_6) 

d^Fo = 8 (7 (V/?i 3 + 4 7 3 /3 2 2 + 12 7 2 /3o/3i/3 2 + 2 7 /3 3 /3 2 + 3 7 /? 1 2 /3 2 + 24/3 2 3 7 /3 « 2 2 + 4/3 1 2 /3 2 2 7 i; 2 2 + (C-7) 

+12 /3i/3 2 2 /3 2 W2 2 + 4 7 2 /3i/3 2 2 w 2 + 12 7 (3 2 (3 2 2 v 2 ~ 6 7 /3 2 /3o/3i 2 «2 - 7 /?i 4 «2 - 3 /3i 3 /3 2 w 2 ^ (C-8) 

/( - 7 2 + 4/3 2 2 w 2 3 -/3i 2 w 2 2 + 2v 27 /3 1 ) . (C-9) 
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